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Abstract 

We discuss the single transverse-spin asymmetry (SSA) to be observed in the D-meson 
production with large transverse-momentum in semi-inclusive deep inelastic scattering, 
ep^ — > eDX. This contribution is embodied as a twist-3 mechanism in the collinear fac- 
torization, which is induced by purely gluonic correlation inside the transversely-polarized 
nucleon, in particular, by the three-gluon correlation effects. The complete formula for the 
corresponding SSA in the leading-order QCD is expressed in terms of the four independent 
gluonic correlation functions and reveals the five independent structures with respect to the 
dependence on the azimuthal angle for the produced D-meson. We present the numerical 
calculations of the SSA formula at the kinematics relevant to a future Electron Ion Collider. 
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1 Introduction 



Charm productions in semi-inclusive deep inelastic scattering (SIDIS), as well as in the pp 
collision, are known to be associated with the twist- 2 gluon distributions in the nucleon, 
since the cc-pair creation through the photon-gluon or gluon-gluon fusion is their driving 
subprocess. Similarly, the twist-3 contributions in the charm productions can be generated 
by the purely gluonic effects inside the nucleon, in particular, the multi-gluon correlations. 
Indeed, the observation of the single transverse-spin asymmetry (SSA) in the open charm 
productions allows us to probe the corresponding twist-3 effects [f, 2, 3, 4, 5]. 

The corresponding SSA arises as a naively T-odd effect in the cross section for the 
scattering of transversely-polarized nucleon off an unpolarized particle, observing a IP- 
meson with momentum Ph in the final state, and this requires, (i) nonzero transverse- 
momentum P^± originating from transverse motion of quark or gluon; (ii) nucleon helicity 
flip in the cut diagrams for the cross section, corresponding to the transverse polarization; 
and (iii) interaction beyond Born level to produce the interfering phase between the LHS 
and the RHS of the cut in those diagrams. In particular, for large P h ± 3> Aq C d, the 
contribution (i) arises perturbatively as the recoil from the hard (unobserved) final-state 
partons, and this leads us to the collinear-factorization framework, so that the other two 
contributions, (ii) and (iii), are generated by the twist-3 mechanism associated with the 
three-gluon correlation functions for the transversely-polarized nucleon [3]. This twist-3 
mechanism may be considered as an extension of the corresponding mechanism for the 
SSA in the pion productions in the SIDIS [10], pp collisions [11, 12], etc., based on the 
quark-gluon correlations in the nucleon, but it has been clarified [3] that a straightforward 
extension [1, 2] to the IP-meson productions leads to missing many terms in the SSA. The 
complete leading-order (LO) QCD formulae for the corresponding SSA in the IP-meson 
productions have been recently derived in [3] for SIDIS and in [5] for pp collisions, based 
on the relevant twist-3 mechanism, i.e., the soft-gluon-pole mechanism, and those results 
revealed, for the first time, the entire nonperturbative gluonic degrees of freedom to induce 
the SSA and the whole structures of the asymmetries with respect to the azimuthal angle 
for the produced IP-meson. Remarkably, those complete twist-3 formulae are related to the 
certain derivative of the twist- 2 cross sections for the IP-meson productions [4, 5], thanks 
to universal structure behind the SSAs in a variety of hard processes [6, 7]. 

The purpose of this paper is to present a numerical estimate of the SSA in the high-P/^ 
IP-meson production in SIDIS, ep^ — > eDX, based on the LO QCD formula of [3]. We use 
the modeling [5] of the twist-3 three-gluon correlation functions in the nucleon, which is 
guided by the SSA for the IP-meson production observed at RHIC [8], and demonstrate the 
influence of the nonperturbative behaviors of gluonic correlations. We calculate the SSA 
at the kinematics relevant to Electron Ion Collider (EIC) [9] . 

2 The LO QCD formula for the D-meson production 

To describe the IP-meson production in SIDIS, e(£) +p(p, S±) — > e(£') + D(P h ) +X, we use, 
as usual, the kinematic variables S ep = (i + p) 2 , q — £ — £', Q 2 = —q 2 , x bj = Q 2 /(2p-q), and 
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Figure 1: Feynman diagrams for partonic subprocess in ep" — y eDX; mirror diagrams also 
contribute. 



Zf = p ■ Ph/{p ■ q)- We work in a frame where the 3-momenta q and p are collinear, both 
moving along the z-axis as q^ = g^ 3 Q and p M = g^Q/ (v^sc&j), an d define qr = Phx/zf and 
the azimuthal angles around the z-axis (ft, $5, and x °f the lepton plane, the spin vector 
Sj_, and the D-meson momentum Pft, respectively [3, 4]. We take into account the masses 
m c and rrih for the charm quark and the D meson. 

In the LO in QCD perturbation theory, the photon-gluon fusion subprocesses of Fig. 1 
drive the SSA for large Ph± D-meson production; in Fig. 1, the above-mentioned contribu- 
tion (i) is provided by the recoil from the hard unobserved c quark and the c quark with the 
momentum p c fragments into the D-meson in the final state. The short bar on the internal 
c-quark line indicates that the pole part is to be taken from the corresponding propagator, to 
produce the interfering phase for the contribution (iii); we note that these pole contributions 
from Fig. 1 would cancel the similar contributions from the corresponding mirror diagrams, 
if the c quark were unobserved in the final state as in the case of the D-meson produc- 
tion. The external curly lines represent the gluons that are generated from the three-gluon 
correlations present inside the transversely-polarized nucleon, (pS j_\A a (0) Ap{rj) A^(£)\pS ±) , 
corresponding to the contribution (ii). The diagrams obtained by the permutation of the 
gluon lines in Fig. 1 also produce the contributions (i)-(iii), but the Bose statistics of the 
gluons in the above matrix element guarantees that we need not consider those diagrams 
separately. Thus, the SSA in the present context can be derived entirely as the contribu- 
tions of soft-gluon-pole (SGP) type [10], leading to k 2 — ki = 0, by evaluating the pole part 
in Fig. 1. 1 The twist-3 nature of those contributions are unraveled by the collinear expan- 
sion, as usual. The expansion produces lots of terms, each of which is not gauge invariant. 
Indeed, many of them vanish or cancel eventually, and the remaining terms can be orga- 
nized into a gauge-invariant form. This can be demonstrated [3] by sophisticated use of the 
Ward identities for the contributions of the diagrams in Fig. 1. The resulting factorization 
formula of the spin-dependent, differential cross section for ep^ — > eDX reads [3, 4] 

d & &° ajnpiselMN /-tt\ v 

[M iWx%si p qA 2 ) k=1 %^ k k 

1 The other types of pole contributions participate in the SSA in other processes, see [10, 12, 13]. 
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(1) 



where [duo] = dxbjdQ 2 dzfdq^d(pdx denotes the differential elements, x = x^/x and z = 
zj/z are the partonic variables associated with the usual momentum fractions x and z, 
respectively, D c (z) denotes the usual twist- 2 fragmentation function for a c-quark to become 
the D-meson, and the quark-flavor index a can, in principle, be c and c, with 5 C — 1 and S 5 = 
— 1, so that the cross section for the .D-meson production ep^ — > eDX can be obtained by a 
simple replacement of the fragmentation function to that for the D meson, D a (z) — > D a (z). 
0(xi,X2) and N(xi,X2) represent a complete set of twist-3 gluonic correlation functions, 
defined through the gauge- invariant lightcone correlation 2 of three field-strength tensors [3], 



M a / ( \(x 1 ,x 2 ) 



( _)V-i,-2, = -^ 3 / ^ / ^e^ 1 e^-* 1 J( P 5 J .|d to Ff , (0)i^(Cn)i!r , (An)|p5 J .; 



--2iM 



N 



0(x 1 ,x 2 )g a ^^' pnS± +0(-x 2 ,x 1 - x 2 )g^e apnS± +0(x 2 - x u -xjg^e 



/?7 apnS± 



■yet J3pnS± 



,(2) 



M 



3 rdA r 
J 2w J 



iC(x2-xi) / 



( P S ± \if bca Fr(0)F2 n ((n)Fr(\n)\pS A 



= 2iM 



N 



N(x 1 ,x 2 )g a(3 e^ pnS ^+N(-x 2 ,x 1 - x 2 )g^e apn ^ +N(x 2 - x u -x^e 



/3y apnS± 



■ya fipnS± 



,(3) 



with the nucleon mass Mn and a lightlike vector n satisfying n 2 = and p-n — 1; d bca and 
jbca are ^ respectively, the symmetric and anti-symmetric structure constants of the color 
SU(3) group, so that 0(x±, x 2 ) and N(xi,x 2 ) are the C-odd and C-even functions satisfying 



0(x 1 , x 2 ) = 0(x 2 , xi) = 0(-Xi, -x 2 ), 
N(xi, x 2 ) = N(x 2 , xi) = -iV(-xi, -x 2 ) 



(4) 



In (1), a em is the fine-structure constant, e c = 2/3 is the electric charge of the c-quark, and 
the summation for the subscript k runs over k = 1,2, 3, 4, 8, 9, with 

Ai = 1 + cosh 2 *0, „4 2 = —2, ^4 3 = — cos(0 — x) sinh 2*0, 

Aa = cos 2(0 — x) sinh 2 -0, ^4 8 = — sin(0 — x) sinh 2-0, Ag = sin 2(0 — x) sinh 2 -0, 

5i = S 2 = S 3 = S 4 = sin($ 5 - x), 5 8 = 5 9 = cos($ 5 - x), (5) 



2 We suppress the gauge-link operators to be inserted in between the field strength tensors. 



4 



where cosh-?/' = 2xbjS ep /Q 2 — 1. The delta function in (1) implies that the lower limits of 
the integrals are given by [1, 3] 



2x bj m 



1 + 



2 

x bjQT 



[1 - x bj ) Q 1 



(6) 



and 



X b j 



X b j 



1 + 



1 + 



z f (l~z f )Q2 



for Zf 1 + A 1 + S >1, 



(7) 



2mf 



1 + A 1 + 



17)2 



for + 4 <1. 



Partonic hard parts Aa l k depend on m c as well as other partonic variables; for the explicit 
formulae of AajL, we refer the readers to Eqs. (71)- (74) in [3]. The participation of the 
"derivative terms" , the terms with the derivatives of the three-gluon correlation functions 
as dO(x,x)/dx, dO(x,0)/dx, dN(x,x)/dx, and dN(x,0)/dx, is characteristic of the con- 
tributions originating from SGPs. Note that, instead of evaluating the SGP contributions 
arising in the diagrams in Fig. 1 as above, those results can be obtained using the "master 
formula" [4], which is schematically given by 



d e A 



a 



[du] 



x(M a / { l ) (x,x)8 a + M a / { l ) (x,x)), 



(8) 



where u^ v = g^ v — p^n v , and H^vixp, q,p c ) represents the partonic hard part for the 2 — > 2 
Born subprocess, expressed by the diagrams in Fig. 1 with the soft (k 2 ~ k± = 0) gluon line 
removed. 3 This reveals that AajL in (1) are related to the twist-2 hard parts 'H lLV {xp, q,p c ), 
similarly as in the SSA in various processes associated with twist-3 quark-gluon correlation 
functions [6, 7]. 

We now reexpress as 



- x = 4>h, 



- <Ps, 



(9) 



in (5), where <ph and <ps represent the azimuthal angles of the hadron plane and the nucleon's 
spin vector S±, respectively, measured from the lepton plane. Then, (1) can be expressed 
as 



d 6 Aa 
[duo] 



= sin(0 h - 4>s) {Fi + T 2 cos <j) h + J" 3 cos 2<f) h ) + cos(4> h - (f> s ) {T± sin <p h + F 5 sin 2<p h ) , 

(10) 



3 Thc color factors of the type, Tr[tH c t a ] = (d bca + if bca )/4, implied by the dia grams in Fig. 1, arc 
contracted with the three-gluon matrix elements, giving rise to the combination of (2) and (3) in (8). 
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with the corresponding structure functions T\,Ti-, . . . , J-5, exhibiting the five independent 
azimuthal dependences similarly as in the twist-3 SSA for ejfi — > enX, generated from the 
quark-gluon correlation functions [6]. Thus, the complete LO QCD formulae (1), (10) for 
the high-P^x -D-meson production in SIDIS are expressed in terms of the four types of 
gluonic functions 0(x, x), 0(x, 0), N(x, x), and N(x, 0) of the relevant momentum fraction 
x, and generate five independent structures about the dependence on the relevant azimuthal 
angles. This is in contrast to the corresponding results in [1], which were expressed by only 
two types of gluonic functions and three independent azimuthal structures. 

The formulae (1), (10) for the single-spin-dependent cross section should be compared 
with the corresponding LO QCD formulae of the twist-2 unpolarized cross section for the 
high-P^x D-meson production in SIDIS, which is generated from the usual unpolarized 
gluon-density distribution G(x), as [3] 

d 6 a unpo1 oLa,e2 A , f 1 dx f 1 dz J q% 1\/ 1\ + m\ \ 



[du] QAn 2 x 2 bj Sl p Q 2 ~ Jx min x Jz min z \Q 2 \ x) \ z) z 2 Q 2 

x £ D a {z)G{x)a u k 



a=c,c 



= + 0% cos 4>h + &z cos 2(f) h , (11) 
where the partonic hard cross sections are given in Eq. (81) of [3] and obey 



! 2q T x v 

Aak = qhT^Y) ° k ' (12) 



as implied by (8). 



3 SSA in the D-meson production at EIC 

We evaluate the SSAs for the D° production, ep^ — > eD°X, based on the QCD factorization 
formula (1). In particular, using (10) and (11), we calculate the asymmetries, 

a?' 2ar 2<t¥ : 2<t? : IoV 

to be observed at a future EIC. Because the partonic hard parts in (1) are common for 
the C-even and C-odd three-gluon functions, we show the contributions to (13) from the 
C-odd correlation functions 0(x,x) and 0(x, 0) in the following calculations. For the first 
estimate presented in this paper, we assume the two types of functional forms of those 
correlations, corresponding to the different small- a; behavior, 

Model 1 : 0(x,x) = O(x,0) = 0.004 x G(x), (14) 
Model 2: Q(x, x) = 0(x, 0) = 0.001 y/xG(x), (15) 
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where the coefficients 0.004 and 0.001 are suggested in [5] by the comparison of the three- 
gluon contribution to the SSA for p^p — > DX with the data observed at RHIC [8]. 4 We use 
the CTEQ6L gluon-density distribution [14] for G(x) of (11), (14) and (15), and KKKS08 
fragmentation function [15] for D a (z) of (1) and (11). For those nonperturbative functions, 
we associate the scale /i 2 = Q 2 + m 2 + zjq 2 ^, and, for simplicity, we assume that the scale 
dependence of 0(x, x) and 0(x, 0) is determined by that of G(x) according to (14) and 
(15). We use m c = 1.5 GeV for the charm quark mass. 

At EIC kinematics with S ep = 1000 GeV 2 , Q = 2.5 GeV, P h± = 2 GeV, and x bj = 0.01, 
the SSAs in (13) are, respectively, shown as a function of Zf in the first five panels in Fig. 2, 
where the contributions due to using the models (14) and (15) in (1) are plotted: the 
solid and dashed curves show the contributions from 0(x, x) and 0(x, 0) of (14) (Model 1), 
respectively, to the asymmetries (13), so that the difference between those two curves reflect 
the difference between the relevant partonic hard parts in (1), i.e., between Aa\ and Act 2 
(Adfc and Aa 4 ); similarly, the dotted and dot-dashed curves show the contributions from 
0(x,x) and 0(x, 0) of (15) (Model 2), respectively, to the asymmetries (13). The behavior 
of aY arising in the denominator in the SSAs (13), as well as of the other coefficients 
cr^3 in the unpolarized cross section (11), is shown as a 1 - = (2irS ep Xbj/ ' zf)crV in the last 
panel in Fig. 2 and demonstrates the D-meson production rate at EIC. 5 We find that the 
contributions to J-'i/cr 1 / are several percent level and significant, while those for the other 
asymmetries are small. With high energy at EIC, the Zf dependence of the asymmetries is 
influenced by the different small-x behaviors between the two models (14) and (15). Similar 
features are observed also in the Ph± dependence of the SSAs (13) and the unpolarized cross 
section (11) for the fixed Zf = 0.3, as shown in Fig. 3. Here the overall behaviors of the 
SSAs obey the 1/Ph± falloff characteristic of the twist-3 effect. A remarkable point revealed 
as a function of Ph± is the growth of the asymmetry T^j (2af ), as well as of the unpolarized 
cross section (11), for decreasing P h± ; in particular, J r 5 /(2crf) could reach a few percent. 
The formulae (1) for the single-spin-dependent cross section tells us that the Ph± as well 
as Zf dependence of the SSAs measured at EIC provides the information to determine the 
x dependence of the gluonic correlation functions 0(x,x), O(x,0), N(x,x), and N(x,0). 
Indeed, the different behaviors of the relevant partonic hard parts in (1) give rise to very 
different shape among the five asymmetries (13) as a function of Zf as well as of Ph±- Some 
asymmetries have maximum and/or minimum, some asymmetries have a node, and some 
asymmetries are monotonic functions. 

The following figures are same as Figs. 2, 3 but for EIC kinematics with higher energies: 
With S ep = 2500 GeV 2 , Q = 4 GeV, and x bj = 0.01, Fig. 4 (Fig. 5) shows the results as a 
function of Zf for Ph± = 2 GeV (as a function of Ph± for Zf = 0.3). With S ep = 5000 GeV 2 , 

4 0{x 1 x) and 0(x, 0) in (14) and (15) are twice as large as those in [5], by which we take into account 
the effect of N(x,x) and N(x,0) contributing constructively as 0(x,x) and 0(x, 0). 

5 Changing the variables as Q 2 — > y (= p ■ q/(p ■ t) = Q 2 /(xbjS ep )), — > Ph±, an d performing the 
integration over \ for a fixed fa, (11) becomes d 5 a unpo1 /(dxbjdydzfdP%j_dfa) — &i cosfa+a^ cos2fa. 
This choice is convenient for comparison with the results in [1]. The results of shown in the last panel in 
Figs. 4 and 5 are reduced compared with the corresponding numerical results in [1]. The reduction mainly 
comes from the use of the quark mass m c — 1.5 GeV corresponding to the pole mass [15], which is larger 
compared with the value m c = 1.3 GeV used in [1]. 
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Figure 2: The individual contributions of (14) and (15) to the SSAs (13) with (10), (1) (the 
first five panels) and the individual coefficients of the unpolarized cross section (11) (the 
last panel), plotted as a function of Zf, for D° production in SIDIS at EIC kinematics with 
S ep = 1000 GeV 2 , Q = 2.5 GeV, P h± = 2 GeV, and x bj = 0.01. 
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Figure 3: The individual contributions of (14) and (15) to the SSAs (13) with (10), (1) (the 
first five panels) and the individual coefficients of the unpolarized cross section (11) (the 
last panel), plotted as a function of Ph±, for D° production in SIDIS at EIC kinematics 
with S ep = 1000 GeV 2 , Q = 2.5 GeV, z f = 0.3, and x bj = 0.01. 
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Figure 6: Same as Fig. 2, but for EIC kinematics with S ep 
P hl _ = 2 GeV, and x bj = 0.005. 
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Figure 7: Same as Fig. 3, but for EIC kinematics with S ep = 5000 GeV 2 , Q = 4 GeV, 
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Q = 4 GeV, and Xbj = 0.005, Fig. 6 (Fig. 7) shows the results as a function of Zf for 
P h ± = 2 GeV (as a function of P h ± for Zf = 0.3). Due to the large values of gluon 
distribution in the small- x region, the cross sections are still sizable for those higher energies, 
and the features found in the SSAs in Figs. 2, 3 are observed also in Figs. 4-7, with some of 
them being even more pronounced. We note that all contributions to the asymmetries from 
the functions 0(x,x) and 0(x, 0), as presented in Figs. 2-7 for the D-meson production, 
change their signs for the Z)-meson production as implied by the factor 5 a in (1), while the 
contributions from N(x,x) and N(x,0) do not. 

The present results with (14), (15) indicates that the derivative terms due to dO(x, x)/dx 
and dO(x, 0)/dx in (1) are the dominant contributions of similar size in the largest asymme- 
try TxjoX , to make its value several percent. On the other hand, for the other asymmetries 
•^2,3,4,5/ (2af), the derivative terms do not give dominant contributions, so that even a non- 
derivative term could give the largest contribution; in particular, the contribution of the 
nonderivative term associated with 0(x, 0)A<5"g in (1) is responsible for the above-mentioned 
growth of J-$/(2aY')i which becomes more pronounced for higher energies. 

For all cases treated in Figs. 2-7, x m [ n in (1) and (11) is given by the formula in the 
first line in (7), so that the value of x m - m becomes the smallest at Zf = 0.5. This prop- 
erty, combined with the fact that the (dominant) derivative-term contributions would give 
J-\/(Ti ~ 1/(1 — rr m in), would explain the existence of the minimum in the solid and dashed 
curves around Zf ~ 0.5 in the first panel in Fig. 2, as discussed in [1]. As demonstrated by 
the dotted and dot-dashed curves in the same figure, however, such behavior is affected by 
the small- x behavior of (15) different from (14), such that the above-mentioned minimum 
around Zf ~ 0.5 could be changed into the maximum. In Figs. 4, 6 with higher energies the 
contributions from the small- x region are more important, and the corresponding minimum 
or maximum is less pronounced. 

4 Conclusions 

In this paper we have discussed the SSAs in SIDIS, ep^ — > eDX, with the D meson having 
large transverse momentum through the photon-gluon fusion mechanism at the twist-3 
level, which is induced by three-gluon correlation inside the nucleon. In particular, we have 
presented a numerical estimate of those SSAs for the first time with the consistent leading- 
order accuracy in QCD, using the corresponding collinear factorization formula. Gauge 
invariance and permutation symmetry among the gluons require the four types of gluonic 
nonperturbative functions and their derivatives to represent the relevant twist-3 mechanism, 
as a consequence of the soft-gluon-pole contributions associated with the 3 — > 2 photon- 
gluon fusion subprocesses or of the master formula with the Born-level (2 — > 2) photon- 
gluon fusion subprocesses. The corresponding SSAs receive the five independent azimuthal 
structures, and our calculation of them, using gluonic nonperturbative functions suggested 
by the RHIC data for p^p — > DX, demonstrates good chance to access multi-gluon effects 
at an Electron Ion Collider, in particular, through the asymmetries with J-\ and T^. The 
similar multi-gluon effects also contribute to the SSAs in ep^ — > enX, as well as in Drell-Yan 
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and direct-photon productions [16]. 
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